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Abstract 

The energy spectrum of two 0-branes for fixed angular momentum in 2+1 
dimensions is calculated by the Rayleigh-Ritz method. The basis function 
used for each angular momentum consists of 80 eigenstates of the harmonic 
oscillator problem on the corresponding space. It is seen that the spectrum 
exhibits a definite linear Regge trajectory behavior. It is argued how this 
behavior supports the picture by which the bound-states of quarks and QCD- 
strings are governed by the quantum mechanics of matrix coordinates. 
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1 Introduction 


The string theoretic description of gauge theories is an old idea EH , still stim¬ 
ulating research works in theoretical physics [4]-[^. Depending on the amount 
of momentum transfer, the hadron-hadron scattering processes have shown two 
different behaviors [^|^. At very large momentum transfers the interactions are 
among the point-like substructures, and qualitative similarities to electron-hadron 
scattering emerge. At high energies and small momentum transfers the Regge 
trajectories are exchanged. The exchanged linear trajectories are the first motiva¬ 
tion for the string picture of strong interaction. However, the fairly good fitting 
between the linear Regge trajectories and the mass of QCD bound-states has not 
been explained yet [^, partially due to the lack of a consistent formulation of 
string theory in 3-1-1 dimensions 

According to string theory, 0-branes are point-like objects to which the strings 
can end [^[l0| . It is known that in a specific regime the dynamics of N 0-branes is 
governed by the matrix quantum mechanics resulting from dimensional reduction 
of U(A^) Yang-Mills theory to 0 -I- 1 dimension 11 . In this regime, the dynamics 
of 0-branes and the strings stretched between them is encoded in the elements of 
matrix coordinates resulted from the dimensional reduction of non-Ableian gauge 
theory. 

By the picture mentioned above, it sounds reasonable that the dynamics of 0- 
branes is used to model the bound-states of quarks and QCD-strings. This picture 
is the main theme of a series of works, and it is shown that the dynamics of 0- 
branes can reproduce some known features and expectations in hadron physics, 
including the potentials between static and fast decaying quarks, and also the 
Regge behavior in the scattering amplitudes 12 13 . The symmetry aspects of the 
picture were studied in 14 . In particular, it is argued that maybe the full featured 
formulation of non-Ableian gauge theories is possible on non-commutative matrix 
spaces 


14 


In the present note the aim is to see whether the 0-brane matrix dynamics can 
reproduce the linear Regge trajectories observed in hadron physics. Early studies 
on spectrum of 0-brane bound-states are reported in 15 -18 . In TopO the study of 
spectrum based on the variational method is presented. In the present work, based 
on the results by TSpO , for the case of two bosonic 0-branes in 2-1-1 dimensions the 
energy eigenvalues are calculated for states with given angular momentum, ranging 
from 0 to 42. The spectrum is calculated by the Rayleigh-Ritz variational method, 
and the basis function for each angular momentum consists of 80 eigenstates of 
the harmonic oscillator problem on the configuration space of the 0-branes. It 
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is seen that apart from two lowest angular momenta, the energy versus angular 
momentum can be htted with straight-line at each level. The spectrum may be 
interpreted as the one for massive 0-branes in 2-1-1 dimensions, or in a Matrix 
theory perspective 21 , as for massless particles in 3-1-1 dimensions but in the 


light-cone frame. Based on the latter way of interpretation, the linear relation 
can be turned as the one between mass-squared and angular momentum, just 
reminiscent the observed one in QCD bound-states. 

Based on the above observation about the spectrum, this may be suggested 
that, the quantum mechanics of matrix coordinates can reconcile string picture 
and QCD in 3-1-1 dimensions. In particular, according to this picture the bound- 
states of quarks and QCD-strings are governed by the quantum mechanics of 


matrix coordinates 12 ■ 14 


The scheme of the rest of this paper is the following. In Sec. 2, the basic notions 
for the 0-brane matrix dynamics are presented, together with a demonstration of 
a bound-state classical solution. In Sec. 3, the quantum theory is developed. The 
eigen-functions of the angular momentum together with the complete solution for 
harmonic oscillator on the 0-branes’ configuration space are presented. This solu¬ 
tion is used to construct the basis function used in the Rayleigh-Ritz variational 
method of Sec. 4. The light-cone interpretation of the results is also presented in 
Sec. 4. 


2 Matrix Dynamics of 0-Branes 

The dynamics of N 0-branes is given by a U(iV) Yang-Mills theory dimensionally 


reduced to 0 -I- 1 dimensions 10 18 , given by (in units h = c= 1) 


L = moTr ^ [A, 

= Dt = dt - i[Ao, ], 


( 1 ) 


with Is as the fundamental string length, and toq = {gsls)~^, with gs as the 
supposedly small string coupling, i.e. mo ^ Y’s are in adjoint representation 
of U(fV) with the usual expansion A = XiaTa, a = The theory is 

invariant under the gauge symmetry 


X ^ X' = UXU\ 

Ao A'o = UAoU^ + iUdtU\ 


( 2 ) 
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where U is an arbitrary time-dependent N x N unitary matrix. Under these 
transformations one can check that: 


DtX -)■ D[X' = U{DtX)U\ 
DtDtX ^ D'^D[X'= U{DtDtX)U^. 


(3) 


For each direction there are variables and it is understood that the extra N'^—N 
degrees of freedom are representing the dynamics of oriented strings stretched 
between N 0-branes. The center-of-mass of 0-branes is represented by the trace 
of the X matrices. 

In the quantum theory the off-diagonal elements of matrices play an essential 
role. In particular, it is shown that in the quantum theory the off-diagonal elements 
cause the interaction between 0-branes. For the case of classically static 0-branes 
it is shown that the fluctuations of the off-diagonal elements develop a linear 


potential, just as the case for QCD-strings stretched between quarks 12 
The canonical momenta are given by: 

Pi = = mo DtXi 

oXi 

by which the Hamiltonian is constructed 

P? 


H = Ti 


_ ^0 \x X U 

2mo 


(4) 


(5) 


As the time-derivative of the dynamical variable Aq is absent, its equation of 
motion introduces a constraint, the so-called Gauss’s law 


Ga . 'y ^ \Xi , i y ^ fabc b Pi c — 1^- 


( 6 ) 


In the present work we take the two dimensional case {d = 2) for a pair of 0- 
branes. It would be quite useful to separate the pure gauge variables from the 
others. For the case of SU(2) theory in 2-1-1 dimensions, following 16 20 we use 
the decomposition 


^ia — {^)ab{A)bj{‘b]^ji ( 7 ) 

in which the matrix is an element of group of SU(2). Accordingly the gauge 
transformations of the variable Xia are captured by through ordinary gauge 
group left multiplications. Parameterizing the SU(2) group elements by the three 
Euler angles, the matrix tk is represented by 

^ = i?,(a)i?,(7)i?,(/I), (8) 
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in which Ra is the rotation matrix about the ath axis. Analogously, the matrix 
r] is an element of the SO(2) group parameterized by the angle 4>, capturing the 


effect of rotation in the two dimensional space. The matrix A takes the form 20 


A = 



(9) 


We mention that the only variable with dimension of length is r. Also, apart 
from pure gauge variables a, /3, and 7 , the two dimensional configuration space 
is spanned by the polar coordinates (r, cj }), and the extra variable 0 appears as an 
internal degree of freedom. 

By the decomposition, the three constraints § take the form [^: 


Gi = sin a cot 'y Pa — sin a esc 7 _P /3 — cos a 
G 2 = cos a cot ^ Pa — cos a esc 7 P /3 + sin a p^ 
G 3 = -Pa 


( 10 ) 


in which p^, and p-y are the conjugate momenta of the pure gauge variables a, 
/3, and 7 . By the constraints ([^, using the explicit forms (10), we have to set: 

Pa = Pp = Pj = 0- (11) 

By imposing the constraints, setting Is = 1 the Hamiltonian takes the form PPol 


^ 2p 




r‘^cos‘^( 20 ) j 8 


( 12 ) 


in which p = m 312 , as the reduced mass appearing in the relative motion of two 
0-branes. It is easy to check that the canonical momentum of 0, p^, is conserved. 
So as expected, the two dimensional angular momentum is a constant of motion. 
The equations of motion by (12) are 

p(f — — 


Pi 


pr3cos2(26») 2 "'^®“ ® 


p{r 6 + 2 f 6 ) - 

1- P<l> 


2 p 2 sin(26l) p 
pr^cos^{29) 2 


sin(26l) cos( 2 d) = 0 


pr"^ cos2(20) 


(13) 


It is easy to check that 0{t) = 0, by which the potential is set to zero, the equations 
for (r, (j)) would come to the form of a free particle in polar coordinate. As an 
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Figure 1: The plots of a numerical solution of (131. The outer curve is representing the 
radial coordinate as a function of the polar angle 4>. The inner one, which is scaled ten 
times to make it visible, is 8{4>). The solution is by the conditions: p = 1/2, Is = 1, 
= 1.42, r(0) = 3, 61(0) = 0.157 rad, r(0) = 0(0) = 0, ./(O) = 0 


illustration that the above equations can develop bound-states, the plots of a 
numerical solution are presented in Fig. 1. In the figure, the outer curve is r((/)) 
as the path of the relative motion of 0 -branes in the polar coordinate setup (r, (j)), 
while the inner curve is a ten times scaled of as the internal degree of freedom 
causing the effective attractive force between 0-branes. Evidently, this solution 
represents an almost circular path for the relative motion of 0 -branes. 


3 Quantum Dynamics 

In passing to quantum theory, the constraints in operator form define the physically 
acceptable states as 

G|V-) = 0 (14) 


By the replacements 


Pa 


. 

^ da' 


PP 


. d . d 


(15) 


one would find, as expected, that the physical wave-functions do not depend on 
the pure gauge degrees of freedom a, /3, and 7 . The Laplacian operator can be 
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constructed using the metric gij 


1 


^ dj) 


(16) 


in which g = det explicitly found to be sin7sin(4^) 20 . So, in the coordi¬ 


nate setup (r, 4^,0), with Q < 9 < tt/A and 0 < (j) < 27r, the Hamiltonian acting 

(17) 


on the wave-function ip{r, 9, (p), takes the form 18 20 

1 


^ = + +^rUm^{29), 


in which 


V?, = 


1 


-d 0 {sm{A9)de) -I- 


dl 


(18) 


^ sin(40)""'’' cos^(20)’ 

Using the scaling ip —>■ the Hamiltonian comes to the form 

^ [^dr (r^dr) + d(V^ - 1V4)) + f rSin2(20), (19) 

By the introduced separation of variables, the two dimensional angular momentum 
, and obviously commutes with the Hamiltonian, [Lz,H] = 0. 


is Lz = -ii^ 0 

So one can construct states with given energy and angular momentum. 

3.1 Angular momentum spectrum 

Here the aim is to find the eigenfunctions and eigenvalues of the operator 

^hyx{o,P) = xyx{9,p) ( 20 ) 


for which we assume as usual 


y^{9,P)=gxi0) 


oirriz ip 


( 21 ) 


with rriz as the quantum number associated to the angular momentum in the two 
dimensional configuration space. Although the pure gauge degrees of freedom have 
been separated out, it is known that a remaining discrete gauge transformation 


would cause that only even integer values are accepted for ruz 20 . In particular. 


the shifts a ^ t: + a and p ^ 2 tt + ip would make equal changes to the original 
variables, namely Xia ^ —Xia, if ixiz has an odd value. So, to construct absolute 
gauge invariant physical states, the quantum number rUz has to be even, setting 


TOz = 2m, m = 0,±l,±2,-' 


( 22 ) 
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Using the change of variable x = cos(40), one has 




(23) 


As the spectrum is invariant under the change m —>■ —m, from now on we take 
TO > 0. Using the replacement gx{x) = (1 + x)‘^^^Qx(x), 


(1 - x'^)Q"{x) + (to - (to + 2)x)Q\x) - ^ 


(24) 


which is known to have Jacobi polynomials of order n = / — to > 0, V^’"^\x), as 

(25) 


solutions 23 . By this the eigenvalue A is found 

A = —16(Z — to/2)(Z — to/ 2 + f), TO < Z = 0,1, • • • 


for the normalized eigenfunction 

yrid,^) = 


21 — m + 1 


{i+cosmr^^vi°j::\cosm) 


^2im(p 


(26) 


2 m+l V y i-m y \ n ^ 

The Jacobi polynomials of our interest satisfy the following recurrence relation, 
which comes mostly helpful when the matrix elements of the Hamiltonian (19) are 
evaluated in the angular momentum basis: 

2(Z + 1)(Z — TO + 1) p(0,m) . X _ 2Z(Z — to) _p(0,m) / x 

(2Z - TO + 1)(2Z - TO + 2) {21- m){2l - to + 1) ’ 

= xVI°jZ\^) ( 27 ) 


(2/-to)(2Z-to + 2)' ' 

3.2 Harmonic oscillator solution 


As we are going to evaluate the spectrum of the Hamiltonian (191 by the variational 
Rayleigh-Ritz method 24 , a set of basis functions is needed, for which we shall 
take those of harmonic oscillator. For a harmonic oscillator with kinetic term as 


in (19) and unit frequency (w = 1), taking 

</) = RE,l,m{r) <(') 

the radial equation would come to the form 


(28) 


1 

2/i 


R 




j™(j- + i) 


RE,l,m J T 2^^ RE,l,m — R RE,l,m (^^) 


in which 


J/" = 4Z - 2TO + 3/2. 


(30) 
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It is known that the above has normalized solutions in terms of the Laguerre 
polynomials 


) — 

with (fc = 0,1, 2, • • •): 


/ 2 fc! 

r(fc + jp + 3/2) 




Ek,i,m = 2k + J/" + 3/2 = 2fc + 4Z - 2 to + 3 


(31) 


(32) 


To calculate the matrix elements of the Hamiltonian (191, the following recurrence 


relations for Laguerre polynomialls would appear mostly useful 23 

4“+')(x)-4“V^(x) = 4“)(x) 

(2A: + a + 1 - x)l[°‘\x) = {k + l)L^“\(a:) + {k + a)L^“\(a;) 

{k + a)L[°‘2^ (x) - (x) = X (a;) (33) 

The following identity for the integral of Laguerre polynomials is known as well [25| 

pOO ^ 

/ 4r" \z)L<i>-^\z)dz = r'l r! x 

^0 

(p + cr)! 


min{( } 

, f7!(k'— a)!(k — a)!(a + t'— k')!((7 + T — k)! 

tT=max{^_“^ } 


(34) 


Of course if max{^_/’} > min{^} the integral is zero. 


k'l 


4 Rayleigh-Ritz Method and Spectrum 

To find the eigenvalues of the Hamiltonian we use the Rayleigh-Ritz vari¬ 
ational method, in which a basis function is needed to approximate the exact 
eigenfunctions. Here we take the basis function to be a collection of eigenstates 
of harmonic oscillator obtained in previous part. As we are interested to find 
eigenvalues with given angular momentum the basis function is taken (recall 
TTiz = 2 m) 

[i’k,i,m,/2ir,9,(j))Y with 

m, , 

!• — 2 ’ ' ' ' ’ 2 ^rnaxi k — 0, • • • , (35) 

in which Umax and determine the level of truncations. By this choice, the 
number of the members of the basis function is equal to (rimax + l)(’2-max + !)■ 
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rriz 

El 

E2 

E3 

E4 

E5 

Ee 

rriz 

El 

E2 

E3 

E4 

E5 

Ee 

0 

2.66 

4.54 

5.95 

7.15 

8.25 

9.09 

22 

13.5 

14.9 

16.5 

18.3 

20.4 

22.8 

2 

4.13 

5.31 

6.22 

7.16 

8.34 

9.79 

24 

14.4 

15.9 

17.6 

19.5 

21.6 

24.1 

4 

5.39 

6.13 

6.89 

7.91 

9.22 

10.9 

26 

15.4 

16.9 

18.7 

20.6 

22.9 

25.4 

6 

6.44 

6.99 

7.83 

8.95 

10.4 

12.1 

28 

16.3 

17.9 

19.7 

21.8 

24.1 

26.7 

8 

7.33 

7.96 

8.89 

10.1 

11.6 

13.4 

30 

17.3 

18.9 

20.8 

22.9 

25.3 

28.0 

10 

8.18 

8.95 

9.97 

11.3 

12.8 

14.7 

32 

18.2 

19.9 

21.9 

24.1 

26.5 

29.2 

12 

9.03 

9.95 

11.1 

12.4 

14.1 

16.1 

34 

19.2 

21.0 

23.0 

25.2 

27.7 

30.5 

14 

9.90 

10.9 

12.2 

13.6 

15.4 

17.4 

36 

20.2 

22.0 

24.0 

26.4 

28.9 

31.8 

16 

10.8 

11.9 

13.3 

14.8 

16.6 

IS.8 

38 

21.1 

23.0 

25.1 

27.5 

30.1 

33.1 

IS 

11.7 

12.9 

14.3 

16.0 

17.9 

20.1 

40 

22.1 

24.0 

26.2 

28.6 

31.3 

34.3 

20 

12.6 

13.9 

15.4 

17.2 

19.1 

21.4 

42 

23.1 

25.0 

27.3 

29.8 

32.5 

35.6 


Table 1: The first six energy eigenvalues for given rUz by the Rayleigh-Ritz method, in 
units For each basis function consists of 80 elements. 


Before to proceed, it would be useful to determine how the spectrum depends 
on the initial parameters Ig and gg (recall mo = l/(3s4)i and g = moj2). By the 


re-scalings 18 


X; 


gy^lgX,, 


(36) 


in the Hamiltonian (|^ one finds that the eigenvalues have the form E = cgl^^ Ig 
with c as dimensionless number (recall we have set h = c = 1 ). 


In calculation of the matrix elements of the Hamiltonian (191 one could avoid 


explicit integrations over r and 9 variables, simply by using the recurrence relations 


(27) and (33), and the integral identity (34). 


For the energy eigenvalues reported in Tab. 1 we have set rimax = ''^'max = 
making 80 elements for the basis function for each niz- 

Apart from two lowest mz’s, the values given in Tab. 1 together with the 
straight-line data fittings are plotted in Fig. 2. The results of the fittings are 


presented in Eq. (37), with the brackets indicating the standard error for each 
given value: 


El = 3.474 [0.059] -I- 0.462 [0.002] m^, 
E 3 = 4.632 [0.020] -I- 0.539 [0.001] mz, 
P 5 = 6.754 [0.038] -|- 0.616 [0.001] m^. 


E 2 = 3.953 [0.031] -|- 0.500 [0.001] m^, 

£'4 = 5.535 [0.027] + 0.579 [0.001] ruz , 

Ee = 8.277 [0.047] -|- 0.654 [0.002] mz- (37) 


By the present standard errors one finds that all the percentage errors are less 
than %2. Further, all the statistical P-values for the straight-line fittings are less 
than 10 “^^, leaving almost no room for the null hypothesis. 
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E , 



Figure 2: The plots of energy eigenvalues versus 
the straight-line fittings given in Eq. (37l. 


according to Tab. 1, together with 


4.1 Light-cone interpretation 

The obtained spectrum may be interpreted as the one for massive 0-branes in 2-1-1 


dimensions, or in a Matrix theory perspective 21 , as for massless particles in 


3-1-1 dimensions but in the light-cone frame. For the latter way of interpretation, 
the relative motion of bound-state constituents (defined by P± = k±i = 0) is 
related to the masses’ constituents and the potential W in transverse directions 


as following 26 


H:=P-=Y: 


‘^Pt 


W 

2P+’ 


in which P+ = and p^’s appear as the masses of constituents but in 

the transverse directions of the light-cone frame. For the case of interest, setting 
rrii = 0 and p^ = p^ one simply has the relation between Hamiltonian eigenvalues 
for relative motion (P_l = 0 ) and the mass squared of bound-state using the key 
relation A4^ = 2P“P+ 21 26 . So in the light-cone frame the linear relation be¬ 


tween the Hamiltonian eigenvalues and the angular momentum turns as the linear 
relation between the mass squared and the angular momentum, just reminiscent 
the observed one in hadron physics. 
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